Solution of inverse problems in elasticity imaging using the adjoint method

We consider the problem of determining the shear modulus of a linearelastic, incompressible medium given boundary data and one component of the displacement field in the entire domain. The problem is derived from applications in quantitative elasticity imaging. We pose the problem as one of minimizing a functional and consider the use of gradient-based algorithms to solve it. In order to calculate the gradient efficiently we develop an algorithm based on the adjoint elasticity operator. The main cost associated with this algorithm is equivalent to solving two forward problems, independent of the number of optimization variables. We present numerical examples that demonstrate the effectiveness of the proposed approach.